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ABSTRACT 

Age derivation techniques for unresolved stellar populations at high redshifts are explored 
using the NUV spectrum of LBDS 53W091 (z = 1.55) and LBDS 53W069 (z = 1.43). 
The photometry and morphology of these galaxies - which are weak radio sources - sug- 
gest they are early-type systems, a feature that makes them ideal test beds for the analysis 
of their ages and metallicities with the use of population synthesis models. In the analysis 
that is based on simple stellar population models, we find a significant degeneracy between 
the derived ages and metallicities both in optical+NIR photometric and NUV spectroscopic 
analyses. This degeneracy is not so strong for LBDS 53W069. However even in this case 
the stellar age cannot be constrained better than to a range roughly encompassing one third 
of the age of the Universe at z = 1.43 (90% confidence level). We have explored several 
independent population synthesis models and consistently found similar results. Broadband 
photometry straddling the rest-frame 4000A break is also subject to a strong age-metallicity 
degeneracy. The use of realistic chemical enrichment assumptions significantly helps in dis- 
entangling the degeneracy. Based on this method, we derive the average stellar age for both 
galaxies around (t*) ~ 3.6 — 3.8 Gyr with better constraints on the youngest possible ages 
(~ 3 Gyr at the 90% confidence level). The comparison with simple stellar population mod- 
els suggest sub-solar metallicities (log Z/Zq = —0.2). A composite model using chemical 
enrichment gives slightly higher metallicities in both galaxies (logZ/Z© = —0.1). Given 
that the stellar component in galaxies forms over times which are larger than a typical chem- 
ical enrichment timescale, we conclude that composite stellar populations must be used in 
all photo-spectroscopic analyses of galaxies. From the observational point of view, the most 
efficient (and feasible) way to set limits on unresolved stellar populations comprises a com- 
bination of Balmer absorption lines along with either low SNR rest frame NUV spectroscopy 
or accurate optical and NIR photometry. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: individual (53W091, 53W069) - galaxies: stellar content 



1 INTRODUCTION 

Estimating the age of the unresolved stellar populations observed 
in galaxies represents a major challenge in our understanding of 
galaxy formation. An ideal observation should allow us to infer the 
star formation history of galaxies from a set of various spectropho- 
tometric observables. The light from a recent burst is dominated 
by the colour of OB stars, whereas the light from old stellar pop- 
ulations come predominantly from G- and K-type giants. Hence, 
as a zeroth order approximation, broadband colours track the stel- 
lar age. However, this stellar clock is not good enough because of 
the effect of metallicity on age estimates so that - within the ob- 
servational uncertainties - the colours of galaxies, whose light is 
dominated by old and metal-poor stars, may be indistinguishable 
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from galaxies mainly composed of young and metal-rich stars. Fur- 
thermore, this degeneracy cannot be simply overcome by a better 
spectral resolution as targeted spectral indices such as the Lick/IDS 
system (Worthey et al. 1994) suffer from an age-metallicity degen- 
eracy similar to broadband photometry (Worthey 1994). 

Accurate spectroscopic dating of stellar populations has been 
attempted over the past years (Stockton, Kellogg & Ridgway 1995; 
Dunlop et al. 1996; Spinrad et al. 1997; Yi et al. 2000) although 
a comprehensive analysis of the effect of metallicity has not been 
considered in detail until recently (Nolan et al. 2003; Ferreras, in 
preparation). 

This paper focuses on LBDS 53W069 (z = 1.43) and 
LBDS 53W091 (z = 1.55), two faint radio galaxies from the 
Leiden-Berkeley Deep Survey (hereafter LBDS; Windhorst et al. 
1984a; 1984b). The search for optical counterparts to faint radio 
emission is a technique which should allow us to spot old stel- 
lar populations at high redshifts (Kron, Koo & Windhorst 1985). 
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However, the light from many of these objects is dominated not 
by starlight but by the active nucleus. An attempt was made to tar- 
get old stellar populations by searching in the LBDS catalog for 
weak radio sources with faint NIR magnitudes (K ^18) and red 
optical— IR colors (R — K > 5). Among the reddest galaxies, 
LBDS 53W069 (K = 18.5; R- K = 6.3) and LBDS 53W091 
(K = 18.7; R — K = 5.8) have been extensively studied (Dunlop 
et al. 1996; Spinrad et al. 1997; Dunlop 1999). Spinrad et al. (1997) 
presented the spectrum of 53W091 using Keck LRIS, which maps 
into its rest-frame NUV (1960 < A/A < 3500). Based on a com- 
parison with their simple stellar populations (SSPs), they found a 
minimum age of 3.5 Gyr, which imposed a significant constraint 
on cosmology. However, due to the substantial difference between 
SSP models (e.g., Yi et al. 2003) and the systematic effects in age 
derivation techniques, the subsequent analyses of Bruzual & Ma- 
gris (1997), Heap et al. (1998) and of Yi et al. (2000) indicated 
significantly younger ages (~ 1 — 2 Gyr). The controversy has con- 
tinued. Nolan et al. (2003) have recently performed a more detailed 
analysis exploring composite stellar populations (i.e. a mixture of 
metallicities) and found ~ 3 Gyr roughly confirming their first age 
estimate. The discrepancy between these two age estimates can be 
translated into a formation redshift. By adopting a ACDM cosmol- 
ogy with Q, m = 0.27 and Ho — 71 km s _1 Mpc -1 (Spergel 
et al. 2003), which is used in this paper hereafter, the literature 
suggests that the stellar component of LBDS 53W091 could have 
been formed at a redshift between zf ~ 2.5 (Yi et al. 2000) and 
z F > 4.7 (Nolan et al. 2003). 

The controversy over the actual age of an allegedly simple 
case such as LBDS 53W091 shows that it is imperative to make 
a robust estimate of the underlying uncertainties. We will ex- 
plore in this paper the uncertainties inherent to any age estimate 
using simple fitting techniques to the observed NUV spectra of 
LBDS 53W069 and 53W091. Given the importance of these es- 
timates to cosmology as well as to galaxy formation, we believe 
the case study presented in this paper is a timely and relevant exer- 
cise which should be borne in mind when extracting ages from the 
integrated properties of stellar populations. 



2 STELLAR DATING USING SSPS 

A first approach to this problem involves the assumption that all 
stars in a given galaxy have the same age and metallicity. These 
so-called simple stellar populations (SSPs) constitute the building 
blocks of population synthesis models, and they give good fits to 
globular cluster data (e.g., Bressan, Chiosi, & Fagotto 1994; Yi et 
al. 2001). Both galaxies targeted in this paper appear morpholog- 
ically to be early-type galaxies with a de Vaucouleurs profile in 
rest-frame B and R bands (Waddington et al. 2002). Their opti- 
cal colours hint at the presence of moderately old stars (Spinrad et 
al. 1997). We have performed a first test comparing the observed 
NUV spectrum with a grid of SSPs over a large range of ages 
(1 <t/Gyr< 4.5) and metallicities (-1 < log Z/Z Q < +0.4). 
The age of the Universe at the redshift of these galaxies (z ~ 1.5) is 
4.35 Gyr with the cosmology adopted in this paper. FigureQshows 
the likelihood contours when performing a \ 2 test, taking into ac- 
count the observed flux and signal to noise ratio in the rest-frame 
range 2000 < A/ A < 3500. The three shaded areas represent 
(from dark to light grey) the 1, 2 and 3<j confidence levels. When 
computing the \ 2 measurements, we normalised both the observed 
and model spectra by their integrated flux in the wavelength range 
specified above. Each panel corresponds to a comparison with dif- 
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Figure 1. \ 2 ma P f° r a comparison of the observed SED of galaxies 
LBDS 53W091 (z = 1.55; left) and LBDS 53W069 (z = 1.43) with sim- 
ple stellar populations from three independent population synthesis models 
(see text for details). In each panel the grey scales correspond to the 1, 2 
and 3a confidence levels from centre to outside. The stars give the position 
of the minimum x' 2 (see TablelT]. 



ferent population synthesis models. From top to bottom: BC02 for 
the latest "pre-STELIB" Bruzual & Chariot (1993) models; YY for 
Yi and Yoon (in preparation) which is an updated version of the Yi, 
Demarque, & Oemler (1997) models, and P2 for Pegase 2 (Fioc & 
Rocca-Volmerange 1997). The stars give the position of the best fit 
for LBDS 53W091 (left) and 53W069 (right), also shown in Ta- 
bleQ along with the best value of the reduced \ 2 ■ The table gives 
the marginalized error bars for the age and metallicity estimates at 
the 90% confidence level. The final number of spectral data points 
in the analysis is 122 (25A resolution) in both galaxies. Figures|2| 
and show the observed SEDs of LBDS 53W091 (Spinrad et al. 
1997) and LBDS 53W069 (Dey et al. in preparation), respectively. 
Three characteristic (la) error bars are shown for reference. The 
inset gives the histogram of signal-to-noise ratios, with the median 
shown by an arrow. The dashed line is the synthetic SED from the 
SSP which corresponds to the best fit from the YY models. The 
bottom panel shows the residuals of the fit as a fraction of the noise 
level. 

One can see that all models give similar results for both galax- 
ies NUV spectral fittings. The figure shows that the confidence lev- 
els are too wide to pin down accurate age estimates; and, thus, the 
wide range of "best" ages and metallicities shown in TableQis to 
be expected. Hence, the most important conclusion one can extract 
from Figure Q is that the NUV spectrum used for the test cannot 
give an accurate age estimate in a simple comparison with SSPs be- 
cause of the age-metallicity degeneracy. This degeneracy is smaller 
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Figure 2. SED of LBDS 53W091 (z = 1.55) show in rest-frame wave- 
lengths (solid line). Three characteristic (1<t) error bars are shown. The 
dashed line is the synthetic SED from the simple stellar population which 
corresponds to the best fit from the YY models (see §2 for details). The 
inset shows the histogram of the observed signal-to-noise ratios. The arrow 
give the position of the median. The bottom panel shows the residuals of 
the fit as a fraction of the noise level. The reduced x 2 that corresponds to 
this fit is shown in TablelTl 



Table 1. Best ages and metallicities (90% confidence level) 



Galaxy 


Model 


Age/Gyr 


log Z/Zq 


X? 


LBDS 53W091 


BC02 


3.07* 


21 
0.8!) 


-0.21° 


14 
0.14 


1.39 




YY™ 


2.87i 


45 
1.11 


-0.15°_ 


34 
0.21 


1.40 




P2 


3.07i 


19 
1.49 


-0.20°_ 


34 
0.14 


1.39 




CSP b 


3.62°_ 


24 
0.53 


-0.09° 


21 
0.32 


1.36 


LBDS 53W069 


BC02 


3.47°_ 


88 
0.G5 


-0.20°_ 


07 
0.0!) 


1.49 




YY 


2.98^ 


43 
0.2G 


-0.14° 


06 
0.20 


1.52 




P2 


3.47° 


84 
0.08 


-0.19° 


07 
0.08 


1.49 




CSP b 


3.84°_ 


31 
0.82 


-0.08°_ 


29 
0.33 


1.43 



a Yi et al. (2000) found roughly 2 Gyr as the best fitting age assuming the 
solar abundance. Here, we used the same models; but we get a larger age 
because lower metallicities are allowed. 

6 Composite models assuming a consistent chemical enrichment history. 
Mass-weighted ages and metallicities are shown. See §4 for details. 



for LBDS 53W069, giving a reasonable constraint on the stellar 
ages. However, notice that at the 90% confidence level, one can- 
not rule out ages which span roughly one third of the age of the 
Universe at the observed redshift. 

Figure|4]further illustrates this point. The inset shows a likeli- 
hood contour similar to those in FigureQfor the YY models. Three 
(age, metallicity) pairs are chosen: #2 corresponds to the best fit, 
#1 and #3 are estimates along the "likelihood ridge". The SEDs 
corresponding to the SSPs for these three points are shown in the 
figure as thin solid, dotted and dashed lines corresponding to points 
1, 2 and 3, respectively. The observed SED of LBDS 53W091 is 
shown as a thick line with three representative la error bars. One 



Figure 3. Same as figure|2]for galaxy LBDS 53W069 (z = 1.43). 



can see the SEDs corresponding to all three points are nearly indis- 
tinguishable even though they span a large range of both age and 
metallicity. 

We have also explored the constraints one could set on a 
comparison with SSPs using the broadband photometric data from 
Waddington et al. (2000). Figure [5] shows likelihood contours in 
age and metallicity when performing a % 2 test with the r — i,i — J, 
J — H and H — K colours. We excluded the Gunn-<? data from the 
analysis because of its large error bars. Through this exercise, we 
also demonstrate the difference between the analysis on the low- 
SNR SED and the one on the broadband photometry with higher 
SNR. The shades (from dark to light grey) correspond to the 1, 
2, and 3<r confidence levels for galaxies LBDS 53W091 (Jeff) and 
53W069 (right), respectively. One can see that the age-metallicity 
degeneracy is stronger than in the case of a comparison with the 
NUV SED. 

Dust was not explored in this comparison. In this sense, these 
age estimates could be considered as upper limits, since the red- 
dening from dust will make a stellar population appear older when 
analyzed by a dustless model. Yi et al. (2000) estimated that a mod- 
est amount of reddening - E(B — V)=0.04 mag - can reduce the age 
by about 0.5 Gyr, reconciling their spectroscopic and photometric 
dating of LBDS 53W091. Furthermore, a small error in the flux cal- 
ibration of the SED will change the overall shape of the spectrum, 
with a significant change in the age and metallicity obtained in this 
way. We want to emphasize here that an accurate spectral dating 
requires a very precise flux calibration. We conclude in this sec- 
tion that a comparison of the available data from LBDS 53W069 
or 53W091 and SSPs over a wide range of ages and metallicities 
cannot give us an accurate value of the age unless an independent 
estimate of the metallicity is obtained. In the next section we ex- 
plore on possible ways to improve the age dating of unresolved 
stellar populations at high redshift. 
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Figure 4. A comparison of the observed Keck LRIS spectrum of LBDS 
53W091 and three simple stellar populations from the YY models (Yi and 
Yoon in preparation). Three representative error bars are shown. The inset 
shows contours of the 1, 2 and 3a confidence levels according to a x 2 test 
as discussed in the paper. Three points in parameter space have been chosen 
to illustrate the wide range of ages and metallicities which are compatible 
with the data. The SEDs corresponding to points 1 , 2 and 3 are shown as a 
thin solid, dotted and dashed line, respectively. 



3 HOW TO IMPROVE THE AGE ESTIMATE 

I. Rest-frame NUV spectroscopy at higher SNR 

Obtaining spectra of high-redshift galaxies requires very long inte- 
gration times on large telescopes. For instance, the rest-frame NUV 
SED of LBDS 53W091 shown in Figure[2]required an effective ex- 
posure time of ~ 20, 000 s using LRIS at the 10 m W. M. Keck tele- 
scope (Spinrad et al. 1997). The spectrum has a median SNR~ 3.5 
per resolution element (25 A). We decided to explore the effect of a 
higher signal to noise ratio on the estimates of the age and metal- 
licity. In the left panels of Figure|6|we generated a synthetic galaxy 
spectrum with the same spectral coverage and resolution as that 
of LBDS 53W091, at SNR= 3 (upper) and SNR= 10 (lower). 
The likelihood map was obtained for a set of 200 realizations of 
the spectrum corresponding to a SSP from the YY models for a 
t = 2 Gyr and Z = Zq (these fiducial values are shown in the 
figure as a dashed line). The shaded contours - from dark to light 
gray - correspond to the 1, 2 and 3a confidence levels. Obviously, 
a more time-demanding spectrum at a higher SNR will reduce the 
width of the "likelihood ridge" although it could not break the age- 
metallicity degeneracy unless very high SNRs (~ 10) are achieved 
- which would imply prohibitively long exposure times even for a 
10 m-class telescope. Hence, for all practical purposes, the accu- 
racy of the age would still hinge on an independent estimate of the 
metallicity. 




Figure 5. Likelihood maps comparing the i, J, H, K photometric data of 
LBDS 53W091 (Jeft) and LBDS 53W069 (right) from Waddington et al. 
(2000) with the predictions for simple stellar populations from four inde- 
pendent modellers (see text for details). The grey shades correspond (from 
dark to light grey) to the 1, 2 and 3<r confidence levels. 



Kauffmann et al. 2003) suggests that a spectral coverage encom- 
passing this break could be a good age indicator. Therefore, we 
decided to test if optical SEDs covering the 4000A break lead to 
more refined age estimates than NUV SEDs. The right panels of 
Figure|6|explore this possibility. We generated 200 synthetic spec- 
tra corresponding to a SSP with t — 2 Gyr and Z = Zq from the 
YY models. These SEDs were computed for SNR= 3 (top) and 
SNR= 10 (bottom). The spectral resolution was assumed to be the 
same as that for the observed rest-frame NUV SED (i.e. AA = 25A 
) and the spectral coverage was chosen to straddle the 4000A break: 
9000 < A/ A < 16500 in the observer's frame, which corresponds 
to 3525 < A/ A < 6450 in the rest-frame. The figure shows that 
one does not achieve a better age constraint at all by using the op- 
tical SED straddling the 4000A break, compared to using the rest- 
frame NUV SED. However, the SNR=10 optical SED yields age 
estimates with markedly better accuracy compared to the SNR=3 
rest-frame NUV SED. For 1-5 Gyr populations, rest-frame optical 
SEDs should take shorter exposures by an order of magnitude to 
achieve the same SNR and but require higher values of SNR. It 
is encouraging to find that high SNR rest-frame optical SEDs can 
be effective tools for deriving ages even at low spectral resolution 
(25A). 



II. Spectroscopy across 4000A break 

The strong age sensitivity of the 4000A break along with its weaker 
dependence with metallicity for young stellar populations (see e.g. 



III. Broadband photometry across 4000A break 

One alternative way of estimating the ages from the flux across the 
4000A break would be to perform broadband photometry. This is 
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Figure 6. Comparison of a synthetic SED from a SSP of the YY models for 
t = 2 Gyr and solar metallicity (represented by the dashed lines). The x 2 
test is performed using a grid of SSPs from the BC03 population synthesis 
model. The contours are shown (from dark to light grey) at the 1, 2, and 3<r 
confidence levels. The synthetic SEDs are generated at the same redshift 
as LBDS 53W091 (i.e. z = 1.552) and for two different signal-to-noise 
ratios, namely SNR=3 (top) - which roughly corresponds to that of the ob- 
served spectrum of LBDS 53W091 (see Figure[2| - and SNR=10 (bottom). 
We have analysed two sets of SEDs corresponding to two different spectral 
regions: near-ultraviolet (NUV, Jeft) with the same spectral range as the ob- 
served one, and optical (OPT, right). The top-left panel also shows the 1, 
2, and 3cr confidence levels of a putative measurement of Balmer absorp- 
tion, namely H/3 ~ 2.8 ± O.lA , (rest-frame), expected for the fiducial age 
(2 Gyr) and solar metallicity assumed for this synthetic SED. See text for 
details. 



much less time-consuming than spectroscopy and it is questionable 
whether one can do better than shown in Figure|5|if a more accurate 
photometry is performed on these galaxies. Figure|7|shows the like- 
lihood map using a synthetic SSP from the YY models for t= 2 Gyr 
and Z — Zq assuming a small photometric error in the colours 
(±0.05 m ). We targeted V - R and R - K colours, which map 
into NUV and optical rest-frame colours, respectively. The shaded 
regions are - from dark to light grey - the 1, 2 and 3<r confidence 
levels. The two upper panels show the analysis when photometry 
is used in the analysis: only V — R colour (top), or when both 
colours (V — R and R — K) are considered (middle panel). This 
approach fares equally well compared to the more time-consuming 
continuum fitting at low SNRs. The age-metallicity degeneracy still 
persists. 



IV. Narrow spectral indices 

After exploring the various approaches discussed above, we are left 
with the option of targeting narrow spectral features. This is a tech- 
nique often used in age estimates of early-type galaxies (e.g. Trager 
et al. 2000; Kuntschner 2000; Bernardi 2003). The key issue is to 
target line indices which have a significantly different dependence 
on age and metallicity. As a simple test, we have checked the like- 
lihood maps one would get from a measurement of a Balmer index 
such as H/3. Balmer indices are very prominent in stellar atmo- 
spheres at T ~ 10, 000 K, which corresponds to main sequence A- 
type stars. Hence, Balmer absorption is especially strong in stellar 
populations over an age range 1 — 3 Gyr, which is ideal for 2 ^ 1.5 
galaxies. Old, metal-poor stars can also contribute significantly to 
Balmer absorption lines (Lee, Yoon, & Lee 2000). However, high- 
redshift galaxies are immune to the complication because they are 
still too young to have developed such stars. The lines in the top- 
left panel of Figure [6] represent the 1, 2 and 3<r confidence levels 
of a simulated measurement of Balmer absorption corresponding to 
H/3 = 2.8 ± O.lA in the rest frame. This is the value to be expected 
for the fiducial model targeted in the analysis throughout §3 (i.e. an 
age of 2 Gyr and solar metallicity). One can see that adding Balmer 
absorption to the analysis helps in constraining the ages. The con- 
tour levels in the bottom panel of Figure^show similar confidence 
levels for a combined measurement of broadband photometry (us- 
ing colours V — R and R — K) along with the same value of H/3 
as above. We assume a higher photometric accuracy ±0.05 mag 
than currently available, i.e. ±0.1 — 0.2 mag (Waddington et al. 
2000). Notice that the use of accurate photometry imposes similar 
constraints on the age and metallicity as in the combined analysis 
of the NUV SED at low signal-to-noise ratios shown in the Figure 
6 top-left panel. Hence, a combination of broadband photometry 
with moderate resolution spectroscopy targeting age-sensitive in- 
dices such as H/3 is the best observational approach to an accurate 
estimate of the age of LBDS 53W091. 

We have chosen H/3 to illustrate the usefulness of Balmer line 
measurements, but higher order Balmer line indices may be pre- 
ferred especially depending on the redshift of the target galaxy. The 
recent study of van Dokkum and Ellis (2003) presents a clear ex- 
ample of such successes. 



4 THE EFFECT OF CHEMICAL ENRICHMENT 

So far, we have based our analysis of stellar ages and metallici- 
ties on simple stellar populations. Age estimates based on SSPs 
rely on the assumption that the stars have the same age and metal- 
licity. Globular clusters are the best candidates for a realistic SSP, 
as inferred by numerous colour-magnitude diagram analyses. Var- 
ious spectrophotometric properties of early-type galaxies hint at a 
very fast star formation process implying a rather small age range 
(e.g. Bower, Lucey & Ellis 1992; Kuntschner 2000). However, the 
radial colour gradients found in elliptical galaxies (e.g. Peletier 
et al. 1990) implies that the stellar component is distributed over 
a significant range of metallicities. Furthermore, low-mass ellipti- 
cals both nearby and at moderate redshifts display rest-frame NUV 
colours (Ferreras & Silk 2000a), Balmer absorption, and [Mg/Fe] 
abundance ratios (Trager et al. 2000), all of which are indicative 
of a more extended star formation history (Ferreras & Silk 2003). 
Hence, comparing early-type galaxies with SSPs may not be such 
a good approximation. 

Furthermore, the morphological analysis of LBDS 53W091 



6 /. Ferreras & S. K. Yi 



o - 



LO 

CD 



o 



LO 

o 



o 



LO 

o 




Age/ Gyr 



Figure 7. The contour plots correspond to an estimate of age and metallicity 
using broadband photometry plus a Balmer index. From top to bottom the 
observational constraints are V — R only; V — R along with R — K or these 
two colours plus H/3 absorption. The photometric uncertainty in this simula- 
tion is ±0.05 mag. H/3 is assumed to be measured with a ±0. lAuncertainty 
in the rest frame of the simulated galaxy (z = 1.55). The stars show the 
position of the best fit. The simulated galaxy has an age of 2 Gyr and solar 
metallicity. 



performed on HST images using WFPC2 and NICMOS in the 
F814W, Fl 10W and F160W passbands is suggestive of a two com- 
ponent system, comprising a compact de Vaucouleurs spheroid 
(r e = 2.7 ± 0.7 kpc) plus an extended exponential disk (h — 
4.2± 1.7 kpc), which only appears in the bluer passband (Wadding- 
ton et al. 2002). This is suggestive of the presence of young stellar 



populations. The contamination in the rest-frame NUV from young 
stars is estimated to be up to 20% (Waddington et al. 2002). 

Therefore, we perform a similar analysis as the one described 
in §2 using a chemical enrichment model which results in a (consis- 
tent) distribution of ages and metallicities. Let us briefly describe 
the parameters which determine the star formation history: 

• The star formation rate (ip) is determined by a power law, with 
a star formation efficiency (Ceff) which is fixed in this analy- 
sis. For simplicity, we decided to use a linear star formation law, 
namely: ip(t) = C'EFFp g (t), where p g is the gas mass. 

• Infall of pre-enriched gas is assumed, at a metallicity Zq/W. 
The infall rate follows a generic "delayed exponential" profile: 
f(t) oc Aiexp(-Ai 2 /2i"|), where At = t - t(z F ), with z F 
being a "formation redshift", and tj is the infall timescale. 

• Gas outflows are parametrized by ^ Bout ^ 1, which 
defines the fraction of the gas returned from stars which is ejected 
from the galaxy. 

The model tracks the stellar, gas and metal components in a sin- 
gle zone. The rest of the details follow standard assumptions about 
stellar evolution and chemical enrichment. We assume a Salpeter 
(1955) initial mass function in the 0.1 < M/Mq < 60 mass range. 
The chemical enrichment model is described in detail elsewhere 
(see e.g. Ferreras & Silk 2000b). With this generic parametrization 
we ran for each galaxy a set of 32 x 32 x 32 star formation histories 
encompassing a wide range of outflow fractions (0 ^ Bout ^ 1); 
formation times (0.4 t(zF') /Gyr ^ 4.2), and infall timescales 
(0.05 ^ r/ /Gyr ^ 2). The upper limit chosen for 77 is motivated 
by the fact that longer infall timescales will result in a significant 
fraction of ongoing star formation at the observed redshift. Hence, 
models with 17 ^, 2 Gyr have a strong NUV component from very 
young stars which is incompatible with the observed SEDs of our 
galaxies. The star formation efficiency was fixed at a high value 
(Ceff = 50 Gyr -1 ) after checking that low star formation effi- 
ciencies were consistently giving higher values for x 2 ■ This value 
of the efficiency is consistent with similar models of star formation 
in elliptical galaxies (Ferreras & Silk 2000b). 

Yi et al. (2000) and Nolan et al. (2003) already considered 
a simple analysis of stellar populations with mixed metallicities, 
although their models did not assume a consistent age-metallicity 
relation obtained by a proper treatment of chemical enrichment. Yi 
et al. (2000) adopted metallicity distributions from instant starburst 
models (and thus no age spread). Nolan et al. (2003) on the other 
hand randomly combined populations from seven metallicity bins 
from 0.01 to 5Zq . In our models, we consistently evolve the metal- 
licity according to the prescriptions described above. It is important 
to note that our star formation history (age-metallicity relations) is 
consistent with the colour-magnitude relations (CMR) of early-type 
galaxies at the present epoch (Bower et al. 1992). In other words, 
our models are not randomly selected but instead calibrated to the 
local CMR information. 

We have performed a x 2 test on these composite models com- 
pared to the observed (SNR~ 3) NUV SEDs of the two galax- 
ies. Figure|8|shows the marginalized likelihood contours as a func- 
tion of the average (mass-weighted) age and metallicity. The solid 
lines give the 1, 2 and 3<r (thick) confidence levels. The best fits, 
uncertainties and reduced chi-squares for these composite models 
are also given in table Q under CSR One can see that the use of a 
chemically consistent model helps in constraining the mean metal- 
licity effectively so that the age estimates can be better determined. 
The 3a confidence level limits the average age between 2.8 and 
4 Gyr, with an average (mass-weighted) metallicity around solar, 
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Figure 8. Likelihood contours comparing the observed NUV SED of 
LBDS 53W091 (top) and LBDS 53W069 (bottom) with a set of composite 
models following a chemical enrichment model (see text for details). The 
contours are at the 1, 2 and 3cr (thick line) confidence levels as a function 
of the (mass-weighted) age and metallicity of the stars. The star gives the 
position of the best fit and the square is the best fit for a comparison with 
simple stellar populations (see §2). The dashed lines are contours of the av- 
erage value of the Balmer H/3 index from 2A , in the upper right corner to 
4A (thick line), in steps of 0.5A . The model predictions fall in the range of 
2 ^H/3/A < 2.5. 



with an uncertainty around ±0.5 dex. Our result is in agreement 
with the randomly mixed metallicity model used in Nolan et al. 
(2003), whose models give a stellar age around 3 Gyr. 

It is important to note that we achieve larger age estimates 
from the NUV spectral fitting analysis when we use composite 
models than when we use SSPs. The star and square in figure [8] 
give the position of the best fit for the analysis using composite and 
simple stellar populations, respectively. The SSP analysis lies at 
the 3a limit of a more consistent approach using a mixture of ages 
and metallicities. Notice that the best fit for a composite model cor- 
responds to an average age around 3.6 Gyr with 0.8Zq (table 0, 



whereas the SSP-based analysis shown in figure Q would require 
~ 2 Gyr if the assume the same metallicity. This result is a con- 
sequence of the mixture of stellar populations with difference ages 
and metallicities, and illustrates the fact that composite stellar pop- 
ulation should be used for these analyses. Furthermore, the best 
chemical enrichment model gives an average age and metallicity 
which would be readily ruled out in a naive approach using SSPs 
(see figure^. 

One effective way to achieve a more accurate age estimate 
lies again in the use of age-sensitive observables such as Balmer 
absorption. The dashed lines in Figure [8] give contours of the hy- 
pothetical H/3 index measurements from 2A - in the upper right 
corner - to 4A (thick) in steps of 0.5A . The chemical enrichment 
model predictions correspond to values of H/3 between 2 and 2. 5 A 
. Hence, Balmer absorption can significantly reduce the uncertain- 
ties in the age estimate and check the validity of the chemical en- 
richment models. Furthermore, we must emphasize here that all 
age estimates based on spectral fittings over a wide range of wave- 
lengths are heavily dependent on a precise calibration of the SED. 
The overall shape of the SED plays a crucial role in all model pre- 
dictions presented in this paper. Therefore, it is very important to 
understand the uncertainties in the flux calibration of the spectra to 
avoid large systematic errors in the analysis. Hence, the analysis of 
Balmer absorption should be considered as a valuable cross-check 
in order to give an accurate answer for the stellar ages. 



5 CONCLUSIONS 

The observed rest-frame NUV SEDs of high redshift weak radio 
galaxies have been claimed to be robust estimators for the ages of 
old stellar populations at high redshifts, which in turn allow us to 
set constraints on the age of the Universe and on cosmological pa- 
rameters (Spinrad et al. 1997). However, in this paper we show that 
the combined effect of age and metallicity results in large error 
bars which are shown to be independent of the population synthe- 
sis model used. This problem persists even if we use the spectral 
energy distribution over a wide range of wavelengths instead of a 
set of broadband filters. Only at SNR=10 or greater can the data 
disentangle the degeneracy. Unfortunately, this corresponds to pro- 
hibitively long exposure times on a 10 m-class telescope. 

A comparison of synthetic SEDs built from SSPs, with noise 
mimicking that of the observed data, shows that signal to noise 
ratios close to those used in the analysis of LBDS 53W091 (i.e. 
SNR~ 3) are not high enough to yield age estimates with apprecia- 
ble precision. Slightly higher SNR will help in reducing the error 
bar, but the age-metallicity degeneracy is still very strong unless 
one can achieve a SNR=10 in rest frame NUV spectroscopy or even 
higher in the optical spectral window. In principle, it may appear 
that the rest-frame NUV would be a desirable window to explore 
for these galaxies since it could pose stronger constraints on the al- 
lowed region of parameter space compared to a rest-frame optical 
SED. However, the weaker stellar continuum and the strong effect 
of dust in this spectral region along with our poorer knowledge of 
stellar emission in the NUV imply that it may be more feasible and 
useful to obtain high precision rest-frame optical photometry strad- 
dling the 4000A break. Old stellar populations are much brighter 
redward of the break, which results in higher SNRs. 

Dating unresolved stellar populations from their spectral en- 
ergy distribution over a wide range of wavelengths is strongly de- 
pendent on the overall shape of the continuum. A small error in the 
flux calibration will distort this shape, thus altering the estimated 
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ages and metallicities. Hence, we want to emphasize that this type 
of studies requires SEDs with a very accurate flux calibration. As 
has been widely known for a decade, the use of Balmer absorption 
lines, in combination with broadband photometry, largely solves 
the problem, and helps us in understanding a possible systematic 
effect derived from uncertainties in the flux calibration of the SED. 
As we discussed in §3, different concerns would lead observers to 
choose different Balmer indices. At z — 1.5 indices such as H/3, 
H7 and H<5 all appear in the NIR spectral window, complicating an 
accurate ground-based spectroscopic measurement. 

Hoping to break the infamous age-metallicity degeneracy 
from a theoretical point of view, we have explored a large set of 
chemically consistent composite population models. Simply be- 
cause no arbitrary metallicity is allowed in such a scheme, this ap- 
proach helped us determine metallicities much better. We found a 
better constraint on the age estimates, giving a range of ages be- 
tween 2.8 and 4 Gyr at the 3<r confidence level for both galax- 
ies, with a metallicity around solar with a ±0.5 dex uncertainty. 
LBDS 53W069 seems to accept models with higher average metal- 
licities. Our results - involving different sets of population syn- 
thesis models and a detailed chemical enrichment scenario - give 
similar results as the analysis of Nolan et al. (2003). Even though 
any model of chemical enrichment introduces further uncertainties 
in the modelling, galaxies should be considered composite mod- 
els as the star formation takes place over times which are always 
longer than the characteristic chemical enrichment timescales. It is 
worth noticing that the best fits for a simple and a composite stel- 
lar population are marginally compatible (see figure|8}. Hence, at 
the expense of adding further uncertainties to the modelling, we 
believe a proper mixture of stellar ages and metallicities should be 
considered in all photo-spectroscopic analyses of galaxies. 
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